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Design Considerations for Radiometer Based on Uncooled 

Microbolometer Detector 

Emest Grimberg - OP GAL 



Abstract 

This paper concludes the large effort sustained by Opgal in designing a radiometer instrument based on the uncooled microbolometar 

A detailed description of the design considerations, temperature drift model and the expected ^accuracy « 
is to enable temperature measurement at a relatively high accuracy for any uncooled m.crobolometer based FLIR, even if the detector 
is not a radiometric one, by using the NUC flag as an extremely low frequency chopper. 

Keywords : FUR, Temperature measurement, Uncooled , Microbolometer .Radiometer, Real time. 
1. Introduction 

1 1 Opgal manufactures a large number of uncooled microbolometer FLIR cameras, covering a large number of 
applications. These cameras contain regular uncooled microbolometer detectors that do not wntam any radiometoc 
shield. A large amount of these cameras comes with a standard lens of 35mm or 50 mm focal length. The mam idea 
is to upgrade these cameras in order to provide radiometric temperature measurement capability under the restriction 
of very minor hardware modifications. 

1 2 The paper presents the studies performed for developing the measurement concept. The paper presents a relatively 
detailed model for the energy exchanged among the microbolometer detector and the view, the optics, and the 
internal camera parts. A mathematical model that describes the detector signal versus target temperature was devised, 
that would result in a polynomial model for temperature measurement The paper also presents the solution 
implemented for reducing the measured temperature drift and the accuracy obtained. There is one major problem not 
Sed in this paper: The two-dimensional modulation transfer function (MTF) influence on the temperature 
measurements has to be compensated by an inverse filter. This subject is only mentioned in this paper but is not 
deeply treated. 



2. Energy exchanged model 

2.1 The video signal at the microbolometer output is generated by the energy exchanged between the detector's 

elements and: 

The view 
« The optics 

■ The internal camera parts. . t , , 

The signal itself is built up as a transient process. Assuming that the readout process is interrupted, the energy exchanged 
takes place until the steady state situation is achieved. The energy exchanged process time constant is equal to the ratio 
between the suspended structure Thermal Capacity and Thermal Conductance. 



2.2 Assuming an unlimited time interval and assuming that the view can be modeled as a Blackbody radiation source 
(emissivity^l) the energy exchanged between the detector and the view is expressed by the following formula: 



P_vieM(ATemd = n A - d f J ] (L{XJem R +ATemP>-L{XJem 1 &yTra 

Equation 1. Energy in joules, exchanged between the detector and the view. 

P view represents the energy [in joules] exchanged between one detector element and the scenery. 
L(X,Temp 0 -h ATemp) represents the view, modeled as an average Blackbody radiation source in Watts/(srm ). 
Temp Q represents the detector's suspended structure temperature in degrees Kelvin. 

* 2 

A det represents the effective area of a single detector's element in m . 
f # represents optics f number. 
X represents the wavelength in meters, 
represents the lower integral limit. 

XJZ represents the upper integral limit. 

Trans _ optics (X) represents optics transmittance 

s(X) represents the emissivity of detector's suspended structure. 

t presents the time in seconds. 

Solving Equation 1 for the time variable, the energy exchanged between one detector's element and the view is described 
by the following expression: 

A 2 

P_ view(ATemp) = - ' - ~ d f * • T(L(/L, Temp^ + ATemp) - L(X,Tem&)) • Trans _optiv{X)z(X) * dX 

4(/#) A _, 

Equation 2. Energy in joules, exchanged between the detector and the view. 

Figure 1 describes the energy exchanged between the detector and the view assuming: 

• f number equals 1 

• high transmittance optics (93%) 

• detector emissivity equals 0.8 

• effective detector area equaling (41*1 0* 6 ) 2 meters 

« a spectral band ranging from 8* 10" 6 to 12.5*10" 6 meters 

2.3The energy P_optics(ATemp) exchanged between the detector and optics assuming similar conditions to those described 
in Paragraphs 2. 1 and 2.2 is given by the following mathematical expression: 

A_2 

P __ optic^ATemp) = 7tr-A_ det- (sin(0)) 2 ](L(X, Temp 0 + ATemp) - Z(A, Temp 0 ))(1 - Trans optic^X))s(X)dX 
Equation 3. Energy in Joules, exchanged between the detector and optics. 
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Figure 1. Energy exchanged between detector and view. 

L{AJemp 0 + Aftap) represents the optics, modeled as an average Blackbody radiation source in Watts/(srm ) and 
emissivity equals (1 - Trans _ optics(Z)) . 

, 1 \ 
0 = arc_ tan( ) 

0 represents half of the planar cone angle related to the light collection. 
6 = 0.464 radian for f#= 1 . 
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Figure 2. Energy exchanged between the detector and optics. 



2 4 The energy P camera(ATemp) exchanged between the detector and the internal camera parts assuming similar 
" conditions to those described in the previous paragraphs 2.1, 2.2, and 2.3 is given by the following mathematical 
expression: 

p_camerc(ATemp) = 7rT'A_d&(\-(sin(0)) 2 ) ]{L(Z,Tem R + ATemp) ~ ^Temp^)) s{X)dX 



x i 



L(A,Temp 0 + ATemp) represents the internal camera parts, modeled as an average Blackbody radiation source in 
Watts/(srm 3 ) and emissivity equals one. 
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Figure 3. Energy exchanged between detector and internal camera parts. 

The energy exchanged between the detector and internal camera parts is the largest one, followed by the energy 
exchanged between the detector and the optics. The energy exchanged between the detector and the view is the smallest 
one Assuming, simultaneously, the same temperature difference between the detector and the view, detector and optics 
and detector against camera internal parts, the amount of energy exchanged between the detector and the view is about 
10%. See figure number 4. . 
2 5 Analysis of the results obtained shows that that the capability to reconstruct the view temperature from the video signal 
' collected is very difficult. The energy exchanged by the detector with the internal camera parts depends on temperature 
distribution of the internal camera parts. This distribution is not constant, due to the fact that the camera's thermal time 
constant is very long (about half an hour) and usually the environmental conditions are not stable for such long periods 
of time A dynamic model based on the temperature distribution of internal camera parts is too complicated and beyond 
the capability to be implemented in the existing hardware. Therefore a different solution has to be adopted. The main 
idea is to use the non-uniformity correction (NUC) update process in order to eliminate the internal camera parts 
influence. The flag inserted in front of the detector during the NUC update process is used as an extremely low- 
frequency chopper. The flag's surface facing the detector was treated and coated in order to behave as close as possible 
as an ideal Blackbody surface. A few very accurate thermistors (±0.025°C) have been installed at the following 
locations: 

2.5.1 One thermistor on the back side of the flag surface 

2.5.2 One thermistor on the external surface of detector's vacuum package 

2.5.3 Two thermistors have been glued to the external optics metal case. 
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Figure 4. View energy exchanged versus total energy exchanged. 



3. Temperature measurement description 

3.1 Let us assume, just for simplicity, that the optics influence is negligible. (Later on we shall review this assumption.) 
Since the NUC update process is periodically performed, we assume that during Process #«, the average detector signal 
As„ (after NUC update and bad pixel replacement (BPR) have been performed) has been calculated. We also assume that 
during the same NUC updating Process number «, the temperature of all the installed thermistors have been acquired. 
Tf n represents the temperature measured on the flag during the a? NUC update process , 

7tf„ represents the temperature measured on the external surface of detector's vacuum package during the number n NUC 
update process, 

To n represents the average temperature measured on the external optics case during the number n NUC update process. 

Let us assume that the average temperature of some arbitrary, relatively large area of interest (the MTF problem will be 
mentioned later on) has to be measured. Let us further assume that the specified area belongs to an ideal Blackbody 
radiation source. 

The first approximation of the measured temperature is expressed by Taylor series: 

were: T _ mes represents the first approximation of the measured temperature. 

Sig represents the average video signal after NUC and BPR of the area of interest mentioned above, at 

some point in time between the n and NUC update processes. 

We define As to be the average detector signal value measured on the flag during the number n NUC updating process 
Therefore by definition: 



T _ mes(As„ ) = Tf„ 



00 



T _ mes = F(S/g - ^ ) + Tf n F(Sig -As n ) = ^a r (Sig -As J 



/=0 

Equation 4. First approximation of the measured temperature 

where: 

F represents the transform that translates the video signal to temperature. 

Equation 4 presents the general idea of using any chopper as a reference source. The chopper's video signal is subtracted 
from the video signal to be measured. In this way a differential measurement is performed. The video value 

Sig — As n is translated to temperature and the offset value Tf n is finally added. This general approach works very well 

in the case that As n and Tf n are almost constant between successive samples. However, as explained above in 

paragraph 2.5 this is not the case. The Error between the real object temperature and the reconstructed one by the first 
approximation formula, was plotted versus the temperature measured on the external surface of detector's vacuum 
package {Td) during a period of more than 10 minutes. See Figure 5. The plot described in Figure 5 shows a typical 
behavior. Figure 6 describes the minimum mean square error linear fit to the Error value. However, the slope of the 
linear best fit is not constant in time. A high slope value means that the camera is far from thermal steady state situation 
and a very low slope means that the camera is in thermal equilibrium. 

3 .2 Each time the NUC update procedure is performed, the error of the measured temperature Tjnes caused by the 

internal camera temperature fluctuations is eliminated. The typical behavior described in Figure 5 and Figure 6 suggests 
that by adding an additional term to the first approximation ( Equation 4) a better solution will be obtained. 

AT wge 

T_mes = F(Sig -As n ) + Tf n + - " (Td - Td n ) 

ATd n 

Equation 5 

where Td n represents the temperature measured on the external surface of detector's vacuum package during the n 
NUC update process. 
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Figure 5. Temperature error versus vacuum package thermistor temperature. 
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Figure 6, Linear best fit to the Error value. 



The general idea expressed in Equation 5 cannot be directly implemented. An approximation of this idea was 
implemented using the values As„,Tf n and Td n> producing Equation 6: 

T _ mes = F(Sig -As n ) + F(As " - 4*3 ~_ j& Z T Ll I {Td - T d„) + Tf„ 

Equation 6. Second Approximation. 

The expression F{As n — As n „ x ) estimates the total flag change in temperature measured by the detector. However this 

value contains the error caused by the internal camera parts change in temperature and the real change of flag's 
temperature. 

3.3 This paragraph deals with F function approximation. A mathematical model that describes the detector's signal after 
NUC and BPR processes was built, in order to estimate the capability to translate the signal to temperature. Figure 7 
describes the output of the mathematical model. This mathematical model was used in order to investigate the accuracy 
that can be obtained using minimum mean square error polynomial approximation to the modeled video signal. For a 
relative small dynamic range of about 60 degrees Celsius a second order polynomial expansion gives an acceptable 
accuracy. For example, a second order polynomial expansion is adequate in order to model the optics contribution to the 
video signal. The third order polynomial approximation covers, with an acceptable error, a range of 300 degrees Celsius, 
while the fourth polynomial approximation spans a range of about 400 degrees Celsius. 




Figure 7. Video output mathematical model 

The minimum mean square error polynomial expansion of order L that samples the data at Q different values, requires 
the solution of L +1 linear set of equations. However, the arithmetic used in order to solve the linear set of equations has 
to be carefully chosen. The minimum number included in the linear set of equations is Q , while the largest number is 

Q 

^L,Si£? L • Tne g ra P n s presented in Figure 8, Figure 9, and Figure 10 have been calculated using high arithmetical 

accuracy. A practical solution might be to divide the temperature range into a number of small regions (including some 
small overlap between the regions), using a lower-order polynomial expansion at each region. Up until this point we 
neglect the optics temperature influence. The next paragraph describes the solution adopted by Opgal in order to deal 
with the optics temperature influence. 
.3 The temperature second approximation described in Equation 6 is based on the existence of function F that 
translates the detector signal to temperature. 



T_mes-Wg-Ai)+ £& ^ 9^ TU (Td-T^+Tfn F(Sig-As n )=% (Sig-AsJ 

Equation 7. 



As was explained in paragraph 2.3 the amount of energy exchanged between the detector and the optics is not negligible. 
Assuming that the optics temperature influence can be expressed by a polynomial expansion of order R and assuming 
that the detector signal translation to temperature can be described by a polynomial expansion of order L, the most 
general form of function F is given by: 
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Figure 8. Second order polynomial approximation. 
The normalized error is shown at left, the temperature error at right. 
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Figure 9. Third order polynomial approximation. 
The normalized error is shown at left, the temperature error at right. 
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Figure 10. Fourth order polynomial approximation. 
The normalized error is shown at left, the temperature error at right 



L R 

F(Sig, To) = ]T a r , - 7b r )5z^ 7 [Temperature in degrees Kelvin] 

Equation 8. 

The minimum mean square error algorithm can calculate the constants a r l assuming that the detector's signal has been 

acquired for different view and optics temperatures. The number of linear equations that have to be simultaneously 
solved is (Z+ \)(R+ 1). Although Opgal uses a 32-bit floating point arithmetic unit, the accuracy is not adequate for such 
a task. Therefore we separate the view temperature from the optics influence. Instead of having one single F function, 
there are several functions F To (Sig) that cover the optics ambient temperature range in small steps of 5 degrees. 



FroiSig) = iX,, • Sig' 

Equation 9 



The optics change in temperature is relatively slow. Consequently, the additional work performed by the digital signal 
processor (DSP) in order to calculate, by linear interpolation, the required a To l constants for the new optics temperature 
is negligible. 

3.5 In order to build a real radiometer, additional variables have to be taken into account, such as the target 
emissivity and ambient temperature influence. However, we would like to mention the influence of the system 
modulation transfer function (MTF) on temperature measurement accuracy. A given object, at a constant temperature, 



will show different temperatures at different distances from the camera. The MTF influence can be compensated, at least 
for some limited spatial frequency range. 




Fourier transform of camera impulse response 

Figure II. Fourier transform of camera impulse response in cycles per mm, including optics detector and 
electronics. 

Figure 1 1 describes a typical two-dimensional Fourier transform of an uncooled microbolometer camera. In order to 
compensate the MTF degradation versus spatial frequency, an inverse filter to the MTF response was designed and 
software implemented. Figure 12 represents the high pass frequency response of the filter designed in order to 
compensate the MTF degradation. 

3.6 The ideas described in this paper have been implemented on an instrument design to detect human beings suffering 
from fever. The instrument is monitoring the face temperature and the areas in the image above some given temperature 
threshold are flickering. The instrument covers a relatively very small range in temperatures between 30 to 42 degrees 
Celsius within an accuracy of ±0.25 degrees Celsius. 
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High Pass Fitter Frequency Response. This filter Is the MTF Inverse 
approximation for a detector that contains 320 by 240 elements. 



Figure 12. MTF inverse filter. 



4 Conclusions 



The basic concept of temperature measurement using a regular uncooled microbolometer detector was described in 
this paper. The additional hardware required for the temperature measurement capability is really limited and 
inexpensive. Opgal expects that the capability to measure the temperature will expand the applications span of the 
regular uncooled microbolometer FL1R cameras. 
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Attached Is a marked version for filing. Please include the changes. 

I did not erase the changes in the claims but I think that they narrow the scope and not needed at this time, 
although they are correct for the final application 

Pis comment «UN OPG 1 draft for filmg_conrected_version,doc>> 

Jonathan Zohn 
Corporate I P Manager 
Elbit Systems group 

This e-mail message has been sent by Elbit Systems Ltd. or ElOp Ltd. and is for the use of the intended 
recipients only. The message may contain information that is client-attorney privileged or confidential. If you 
are not the intended recipient you are hereby notified that any use, distribution or copying of this 
communication is strictly prohibited, and you are requested to delete the e-mail and any attachments and 
notify the sender immediately. 



This e-mail message has been sent by Elop Electro-Optics Industries Ltd. and is For the use of the 
intended recipients only. The message may contain Information that Is privileged or confidentiaLIf 
you are not the intended recipient you are hereby notified that any use, distribution or copying of 
this communication is strictly prohibited, and you are requested to delete the e-mail and any 
attachments and notify the sender immediately. 



07/08/2003 



Message 




Geoff Mel nick 



From: Jonathan Zohn [ELOP] [el14012@elopxo.il] 
Sent: 19:33 2003 oouik 07 -wnn oi« 
To: 'Geoff Melnick' 
Subject: RE: OPG 1 

I really had something like that in mind, 
claim 3: a flag will be better than the flag 

OK TO FILE IN THE IPO 

in the name of OPGAL Ltd. 

Jonathan Zohn 

I P Manager 



Page 1 of 2 




This e-mail message has been sent by Eibit Systems Ltd. or ElOp Ltd. and is for the use of the intended 
recipients only. The message may contain information that is client-attorney privileged or confidential. If you 
are not the intended recipient you are hereby notified that any use, distribution or copying of this 
communication is strictly prohibited, and you are requested to delete the e-mail and any attachments and 
notify the sender immediately. 

Original Message 

From; Geoff Melnick [maifto:geoff@ipatent.co.if] 
Sent: Thursday / August 07, 2003 7:23 PM 
To: Jonathan Zohn [ELOP] 
Subject: RE: OPG l 

Hi Jonathan, 

Attached please find my revisions to the claims. 

What I have done is removed the additions from claim 1 and written them out as separate dependent 
claims. That way you get the best of both worlds. It is clearly written out but doesn't limit. 



Geoffrey. L. Melnick 
Patent Attorney 
EHRLICH & PARTNERS • 
O. E.Ehrlich (1995) Ltd. 
Ayaloa Tower, 15 th Floor, 

1 1 Meaacbem Begin Street, 52 521 Ramat Gan, Israel 

Tel.: 03-6127676 Fax.: 03-6127575 Mobile 054 962625 Home 08 936 2304 
My address: geoEf@ipatentco.ti 
Our Website: wvAv.ipateriLco.il 

Alliance: Granot, Strauss, Adar & Co. - Website: www.granot-law.co.il 

This e-mail message is for the sole use of the intended recipients) and may contain confidential and privileged information. Any 
unauthorized review, use, disclosure or distribution is prohibited. If you are not the intended recipient, please contact the sender by 
reply e-mail and destroy all copies of the origind message. 
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intended recipients only. The message may contain information that is privileged or confidential^ 
you are not the intended recipient you are hereby notified that any use, distribution or copying of 
this communication is strictly prohibited, and you are requested to delete the e-mail and any 
attachments and notify the sender immediately. 
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